Probing two-fluid hydrodynamics in a trapped Fermi superfluid at unitarity 
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We develop a variational approach to calculate the density response function at finite temperatures 
of the lowest-lying two-fluid modes in a trapped two-component Fermi superfluid close to a Feshbach 
resonance. The out-of-phase oscillations, which are the analogue in trapped gases of second sound 
in uniform superfluids, have so far not been observed in cold-atom experiments. At unitarity, we 
show that these modes are observable at finite temperatures via two-photon Bragg scattering, whose 
spectrum is related to the imaginary part of density response function. This provides direct evidence 
for superfiuidity and a promising way to test microscopic results for thermodynamics at unitarity. 
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There is growing experimental evidence [H, [3i 01 that 
trapped Fermi gases close to a Feshbach resonance are 
in the collisional hydrodynamic regime. It has been ar- 
gued that the finite temperature dynamics in this region 
should be described by the Landau two-fluid hydrody- 
namic equations Landau's two- fluid hydrodynam- 
ics 0, [gII describes the coupled dynamics of the superfluid 
and normal fluid components when collisions produce a 
state of local thermodynamic equilibrium. A Feshbach 
resonance can be used to increase the magnitude of the 
s-wave scattering length between Fermi atoms in two dif- 
ferent hyperfine states. Thus, close to the unitarity limit 
where the scattering length diverges, the Landau two- 
fluid description of finite temperature collisional dynam- 
ics of trapped Fermi superfluids is expected to be valid. 
In a recent paper Q , we have used a variational approach 
to calculate the two- fluid dipole and breathing modes. In 
this Letter, we extend this variational theory to calculate 
the density response function given by Landau's two-fluid 
hydrodynamics in a trapped Fermi gas at unitarity. We 
propose that two-photon Bragg scattering [8|l (which has 
been used to study the coUisionless region [1, 13, lH) is 
an ideal tool to probe Landau two-fiuid behaviour. 

The most spectacular prediction of Landau's two-fluid 
hydrodynamics is the existence of second-sound, an out- 
of-phase oscillation of the normal and superfluid compo- 
nents. Due to weak interactions, this region has never 
been probed in Bose-condensed gases. For a strongly in- 
teracting Fermi gas near unitarity, one can show that 
the recent studies of breat hing modes have probed the 
in-phase two-fluid mode [J, Il2l. [ist since only this mode 
is excited by modulating the harmonic confining poten- 
tial. We have shown in Ref. 0] that the frequency of 
this in-phase mode is temperature independent at uni- 
tarity. This result explains why these experiments arc 
well explained by the oscillations of a pure superfluid at 
T = [3, m, The out-of-phase modes predicted 
by two-fiuid hydrodynamics in trapped Fermi gases at 



finite temperatures are much more interesting, being the 
analogue of second sound. 

The Bragg scattering cross section is proportional to 
the imaginary part Imxpp(q, f^) of the density response 
function, where q and Q are the wave vector and fre- 
quency difference between two laser beams. In experi- 
ments done so far on trapped Bose gases, Bragg scat- 
tering has been used both in the large q (= |q|) region 
{qRtf ^ 27r, Rtf being the smallest Thomas- Fermi 
radius of the condensate) where the phonon excitations 
characteristic of uniform systems are probed llj , and in 
the small q region such that qExF ^ 0{l) In the lat- 
ter region, the wavelength \ = 2'k j q of the Bragg beams 
is comparable to the size of the condensate, and therefore 
the lowest energy normal modes can be excited. 

The variational approach we develop (see also Ref. 
is well suited for dealing with the two-fluid hydrody- 
namics at finite temperatures in a spatially non-uniform 
trapped Fermi gas. We concentrate on the BCS-BEC 
crossover at unitarity since this is perhaps the most in- 
teresting region. However, our formalism is applicable 
away from unitarity as long as the collisions are strong 
enough to produce local equilibrium. Our new results 
show that the density response spectrum exhibits well 
defined resonances corresponding to the out-of-phase os- 
cillations at temperatures of the order O.STc in the su- 
perfiuid phase. The frequencies and spectral weights of 
the out-of-phase modes are very dependent on the tem- 
perature. The weight vanishes both at T = and above 
Tc- The in-phase hydrodynamic modes frequencies are 
independent of temperature 7] but the spectral weights 
change with T. 

Out-of-phase mode frequencies. — Before discussing 
the density response function, we need to review the vari- 
ational theory [1, 0] for the solutions of the Landau two- 
fiuid equations. This is based on an action that 
describes fluctuations in the displacement fields Us,m„ 
deflned by v.; = lis and v„ = iin, where Vj and v„ are 



2 



the superfluid and normal fluid velocities, respectively. 
The stationary point of the action gives the solutions of 
the linearized Landau two-fluid differential equations and 
hence, the two-fluid hydrodynamic mode frequencies. 

Following Refs. [1, 0|, we make a physically moti- 
vated ansatz for the displacement fields for the low en- 
ergy hydrodynamic modes in trapped gases based on the 
known exact solutions for a pure superfiuid (T = 0) 
and a normal fiuid (T > Tc). In these two limits, the 
exact solutions are of the same form. For the dipole 
mode, characterized by displacements of the centre-of- 
masses of the normal and superfluid components along 
one of the axes of the harmonic trap (we choose the 



z axis), we use Us,Ti(r,t) = J2u 



as,„(w)z, where z 



is the unit vector along the z axis. For the breathing 
mode in an isotropic trap, we use the scaling ansatz 
Ws_„(r, i) = ^'^^ both these modes, 

the action can be written in terms of the variational pa- 
rameters as and a„ as 



^(2) 



where we have deflned the 2x2 matrix A as 

,2 
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and the superscript "T" denotes the transpose of the two 
component spinor. The superfluid and normal fluid mass 
moments, Alg and M„, as well as the "spring constants" 
fcs, fc„, and ksm depend on the choice of ansatz for dif- 
ferent normal modes [1, 0]. These are given below for 
the dipole and breathing modes. The variational solu- 
tions (denoted by the superscript "L") of the Landau 
two-fluid equations are given by 



dS, 



(2) 



da* 



ds: 



(2) 



dal 



0, 



(3) 



which leads to A(af',a^)-^ — 0. The frequencies of the 
in-phase (wi) and out-of-phase (^2) modes are thus given 
by the solution of detA(a;) = 0. As shown in Ref. [7!|, for 
both the dipole and breathing modes at unitarity, the 
in-phase and out-of-phase modes are described by the 
displacements as = and Msas + Af„a„ = 0, respec- 
tively. 

For the dipole mode, one finds [7] AI^ = J dr pso(r), 
^Jdr p„o(r), kg = M^Lol - kg and k^ = (M^ - 
Mf )wg -I- fcf , where we have defined 



k 



D 



dr 



dp J ^ dz dz 



(4) 



Here, psQ and pnQ ^ Pa — Pso are the equilibrium su- 
perfluid and normal fluid densities, respectively. Using 
these results, the in-phase dipole mode is just the gener- 
alized Kohn mode with frequency ujid = equal to the 



trap frequency at all temperatures. The frequency of the 
out-of-phase dipole oscillation is 
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where Mg = MgM^/{Mg + Af,f ) is the dipole mode 
reduced mass. 

At unitarity, the universal nature of the thermodynam- 
ics yjl] leads to great simplifications in the expression for 
the breathing mode frequencies. For an isotropic trap, 
we find 0], Mf = / dr pso{v)r^ , M,f = / dr p„o(r)r2. 



kg = AMgiol - kg and k^ 
where we have defined 



4(M„^ - Mg)ujl -f ki 



dr 



dp. 
'dp 
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The frequency of the in-phase mode at unitarity is inde- 
pendent of temperature, given by ujib — 2wo 'iS!]. The 
frequency of the out-of-phase mode is 
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(7) 



with the breathing mode reduced mass A'lg = 
MfMg/{Mg + Mg ). The frequencies of both this out- 
of-phase breathing mode and the dipole mode [Eq. ([7])] 
are strongly dependent on temperature (see Figs. 3 and 
4 of Ref. 0). 

Density response function. — The two-photon Bragg 
scattering spectrum is well known to be directly related 
to the imaginary part of the density response function 
Xpp (q, fi) [111 ll9(| . This in turn is related to the antisym- 
metric part of the dynamic structure factor S'(q, f2), 



Imxpp(q, ^) 
2TT7n^N 



-5(q,-n)] = 5A(q,f]), (8) 



where iV is the number of Fermi atoms. This formalism is 
valid in both the coUisionless and hydrodynamic regions. 
In this Letter, we concentrate on the contribution to Sa 
from the hydrodynamic dipole and breathing modes in 
trapped gases. 

To find the density response function Xpp 
we consider an external potential of Bragg beams 
which excites modes of frequency VL and wavevector q: 
l^crt(r,i) = V^j^ae*^'! '""^*'. Assuming this is a weak per- 
turbation, there will be a density response of the form 



(9) 



which defines the linear density response function Xpp- 
Our task is to find the structure of Xpp (Qi ^) imposed by 
the Landau two-fiuid equations in the presence of both 
V^crt and the trap potential. 

The Bragg perturbing potential generates an extra 



term S, 



(2) 
pert 



J drdt[V^ort(r, i)5jO+(r,t)-)-h.c.] in the ac- 
tion, with 5p{Y,t) = -V-[pso{r)us{r,t) + pno{r)un{r,t)] 
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being the density fluctuation operator Thus, 

41 = K,,o [F:{q)a:{n) + F„*(q)a:(r!)] + h.c, (10) 
with the weighting factors 

F,^„(q) = -iq^J dre*'i >,o.„o(r) 

F,U^) ^ -ij dre'i-lq-Op.o.noW (11) 

for the dipole (D) and isotropic breathing modes (B). The 
variational solutions of the Landau two-fluid equations 
are now given by Eq. ([3]), using the action S'q^'' + S'pg'j.^.. 
It is straightforward to show that this gives (af', a^)^ = 
— A~-^ [F*(q), F*(q)]"^ Vq_n for the variational parame- 
ters describing the two-fluid modes. Hence, the density 
response of the Landau two-fluid equations, ^PL(q, ^) — 
-[Fs{q)as{Vl) + Fn{q)an{Vl)\, is given by Eq. P but 
now we have an explicit expression for the density re- 
sponse function in the Landau two-fluid region, 



X^p(q,f^) - ^^^^^:[A-i(f^)]„^F^, 



(12) 



where \A'^^{^)]ai3 (a,/? = s,n) denotes the a(3 element 
of the inverse of the matrix A{U) deflned in Eq. 

Using Eq. ([2]) in Eq. (fT2|l . the density fluctuation spec- 
trum given by Eq. ([8]) becomes (for > 0), 



5f(q,0)=7q 



2uJi 



^2(q) 

2^2 



where 



1 



7q 
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,(13) 



(14) 



and the relative weight of the in-phase mode wi is given 
by 



^i(q) 



\Fs{q) + FnjqW Mr 
\Fs{q)\^ Mn + \Fn{q)f Ms 



(15) 



The functions i^„,s(q) are defined in Eq. pT|) . The 
out-of-phase mode uj2 in Eq. (fT^]) has relative weight 
Z2{q) = 1 - Zi(q). Equations ^ and ^ are valid 
for both the dipole (with Fs,n = F^„ and Ms,n = M^J 
and breathing {Fs,n = F^^i and Ms^n = Mf^^) modes. 
Equation pS]). describing the density response function 
predicted by Landau two-fluid hydrodynamics, is the key 
result of this Letter. It is the natural generalization of 
the well-known expression [2^ for a uniform superfluid, 
to be discussed shortly. 

Numerical results for (q, il) for the dipole and 
breathing modes for several different temperatures are 
shown in Fig. [TJ We use the thermodynamic functions 
and superfluid density for a uniform superfluid given in 
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FIG. 1: (Color online) Temperature dependence of the den- 
sity response spectrum as a function of the frequency Q. at 
unitarity due to (a) dipole and (b) breathing modes, for an 
isotropic trap with frequency ujq. The momentum transferred 
q — 2-k/ Rtf is given in terms of the Thomas-Fermi radius of 
a ideal trapped Fermi gas. From bottom to top, the temper- 
ature T increases from 0.04Tf to 0.28rF with a step 0.04rp, 
where Tf is the Fermi temperature of an ideal trapped Fermi 
gas. The top curve is above the superfluid transition temper- 
ature Tc ~ Q.27Tf. 




FIG. 2: (Color online) The spectral function SA(q,<i') with 
increasing Bragg momentum q, for T = 0.7Tc. 



Ref. 01 and use a local density calculation approxima- 
tion to calculate these quantities in a trapped gas. The 
results are quite dependent on the momentum transfer 
q = 27r/A for the Bragg pulse. This is shown in Fig. [21 
which compares S^{q,Vl) for several different values of 
q, at the temperature T — O.TTc. In Figs. [1] and [21 the 
width of the resonances is taken to be 0.02ti;o. In experi- 
ments, the width will go as 27r/Ai, where At is the Bragg 
pulse duration (see, for example, Ref. _19]). 

As T^O, Mr^Mn, and the weights of both the in- 
phase dipole and breathing modes goes to unity, Zi 
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FIG. 3: Momentum dependence of (a) the /-sum weight 7q (in 
units of /Am) and (b) weight Z2 of the lowest energy out- 
of-phase dipole and breathing modes at T = O.ISTf — O.TTc. 



F2(q)A//jF2(q)Af„ = l. Thus, S'f (q,w) only has reso- 
nances at the frequencies of the in-phase modes at T = 0. 
As the temperature increases, however, the out-of-phase 
modes acquire finite weight and S^{q,il) gains addi- 
tional resonances at these frequencies. As T — s- Tc, 
Mr ^ Ms, and Zi -> F^{q)Mr/ F^^{q)Ms = 1. As ex- 
pected, above Tc, only the in-phase modes are present. 
In understanding ^^(q, f^) in both uniform and 



trapped gases, the well-known /-sum rule [8|, |2l| is very 



useful. Quite generally, 5'yi(q, f^) satisfies the condition 



(16) 



The spectral density Sa has contributions from all hy- 
drodynamic modes, as well from high energy excitations 
not described by the two-fluid equations. For a trapped 
gas, the two-fluid density response in Eq. (|13p gives 



/•OO 

Jo 



(17) 



which shows that 7q in Eq. (fT4|) describes the relative 
contribution of a hydrodynamic mode to the /-sum rule 
for a given value of q. 

In Fig. [21 we plot the dependence on q of the weight 
factor 7q in Eqs. p3)) and ([15]) of the dipole and breath- 
ing hydrodynamic modes as well as the weight ^2(q) of 
the out-of-phase modes. These results are for T = 0.7Tc 
and are quite dependent on temperature (as can be seen 
in Fig. [T|). The optimal momentum to observe the out- 
of-phase mode is seen to be about 2-k/Rtf, where the 
overall weight 7qZ2(q) reaches a maximum. As shown in 
Fig. [Sl^a) , 7q is significantly reduced from the maximum 
value of /Am when q ^ 2tt/Rtf- This means that 
at larger values of g, the spectral weight in ^^(q, fi) is 
shifting to other low-energy hydrodynamic modes (e.g., 
quadrupole) as well as high-energy excitations not de- 
scribed by two-fluid hydrodynamics. The successful ob- 
servation of the out-of-phase modes using Bragg scatter- 
ing will require trying to find the optimal values of q and 
T. 




FIG. 4: (Color online) Temperature dependence of the den- 
sity response spectrum in a uniform gas as a function of the 
frequency SI. The offset spectra go from T — OAT^ to T — Tc 
in steps of O.lTc. The spectral weight of second sound (lower 
branch) is very small below OATc- The spectral weight of first 
sound is everywhere much larger than second sound and has 
been reduced by a factor of 5 in this plot. 



Using a plane- wave ansatz for the displacement fields, 
one can verify that Eq. ([T^) gives the usual two-fiuid 
density response function for a uniform superfluid [20| . 
In Fig. [H we show S'^(q, SI) for a uniform gas at uni- 
tarity, based on the same thermodynamic functions 
used to compute the density fluctuation spectrum for a 
trapped gas (using a local density approximation) shown 
in Figs. [1][31 For a uniform gas, the poles of S^{q,n) 
correspond to first and second sound {ui — ciq), the ana- 
logue of the in-phase and out-of-phase modes in a trap 
which we have been discussing. In a uniform superfiuid, 
the expression for S'^(q, fi) describing first and second 
sound modes saturates the /-sum rule in Eq. (|16p , namely 



/•OO 

/ df7 llS'f (q,0) = 
JO 



q^ /Am. 



(18) 



However, the nature of both the in-phase and out-of- 
phase modes is quite different in traps and there is no sim- 
ple comparison. In contrast to the results for a trapped 
gas shown in Fig. [TJ the second sound frequency goes 
to zero at Tc in a uniform gas. The lowest temperature 
shown in Fig.[5|is 0.4Tc, since the weight of second sound 
is vanishingly small at lower temperatures. However, we 
note that the second sound velocity C2 increases to the 
T = value c\/\fZ ~ 0.21i;_f, where C\ ~ O.STuf is equiv- 
alent to the Bogoliubov phonon velocity at T = (see, 
for instance, Ref. (2^). 

Conclusions. — We have proposed that two-photon 
Bragg scattering is an ideal tool to probe the Landau 
two-fluid superfluid hydrodynamics in a trapped strongly 
interacting Fermi gas near unitarity. To this end, we have 
presented explicit predictions for the density response 
spectrum which can be measured with available exper- 
imental techniques [^, [l^, 0| . Previous studies for uni- 
form 17 , 22 and trapped [231 gases did not calculate the 
spectral weights of the hydrodynamic modes in the dy- 
namic structure factor. The observation of second sound- 
like out-of-phase modes at finite temperatures is one of 



the challenges in current research on ultracold atomic 
gases. 
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